Hyperfine interaction and electron-spin decoherence in graphene and carbon 

nanotube quantum dots 
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We analytically calculate the nuclear-spin interactions of a single electron confined to a carbon 
nanotube or graphene quantum dot. While the conduction-band states in graphene are p-type, 
the accordant states in a carbon nanotube are sp-hybridized due to curvature. This leads to an 
interesting interplay between isotropic and anisotropic hyperfine interactions. By using only ana- 
lytical methods, we are able to show how the interaction strength depends on important physical 
parameters, such as curvature and isotope abundances. We show that for the investigated carbon 
structures, the 13 C hyperfine coupling strength is less than 1 /xeV, and that the associated electron- 
spin decoherence time can be expected to be several tens of microseconds or longer, depending on 
the abundance of spin-carrying 13 C nuclei. Furthermore, we find that the hyperfine-induced Knight 
shift is highly anisotropic, both in graphene and in nanotubes of arbitrary chirality. 

PACS numbers: 03.65.Yz, 72.25.Rb, 73. 21. La, 31.30.Gs, 61.48.De 



I. INTRODUCTION 



The spin of an electron confined to a semiconductor 
quantum dot is a prime candidate for quantum infor- 
mation processing devices, with potential applications in 
spintronics 1 ^^ and quantum computational^ Schemes 
for quantum computation rely on a sufficiently long life- 
time of initialized spin states. One major problem is 
that the electron is not isolated from its environment 
but interacts with the nuclei in the semiconductor it has 
been confined to. This interaction leads to relaxation 
of excited spin states, as well as to the decay of spin- 
state superpositions (decoherence). For most quantum 
dots at low temperatures, the main source of spin de- 
coherence is the electron's interaction with a fluctuating 
magnetic field created by the randomly precessing nu- 
clear spins jl i 8 ! 9 ! 10 ! 11 ! 12 ! 13 For III-V semiconductor quan- 
tum dots, with their high abundance of spin-carrying nu- 
clear isotopes, the associated decoherence times are quite 
short, usually of order nanoseconds ^ 9 ' 10 : 14 if no manip- 
ulations on the nuclear system are performed. 

One way of overcoming the problem of short deco- 
herence times is to build quantum dots from semicon- 
ductors with lower abundances of spin-carrying isotopes, 
potentially resulting in weaker nuclear-spin interactions. 
Carbon structures naturally consist of 99% 12 C with nu- 
clear spin and only of 1% 13 C with nuclear spin |, 
and are therefore promising materials for building quan- 
tum dots featuring long spin decoherence times. This 
extraordinary property of carbon materials, as well as 
their comparatively weak spin-orbit interactions, has led 
to proposals of fabricating quantum dots in graphene 
ribbons with armchair boundaries 15 and in carbon nan- 
otubes (CNTs)J£ 

Although it was not yet possible to experi- 
mentally realize a single-electron quantum dot 
in graphene, Coulomb-blockade measurements 
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FIG. 1: Sketch of the systems under consideration through- 
out this work: (a) a single electron confined to a quantum dot 
defined by two barriers on a graphene ribbon; (b) a single- 
electron quantum dot defined in a carbon nanotube. Exter- 
nally applied magnetic fields are indicated by the red arrows. 



in gated graphene quantum dots have already 
been carried out successfully! 17 ' 18 Experiments 
on single electrons in CNTs have advanced even 
morc _ 19.20,21.22,23.24,25.26, 27.28.29.30,31.32,33 Very re- 
cently, first measurements on nuclear-spin interactions 
and electron-spin dynamics in CNTs have been 
performed ) 34 ' 35 reporting an unexpectedly strong hyper- 
fine interaction of order A ~ 100 fieV in 13 C enriched 
CNTs (we will comment on the discrepancy between 
this result and our calculations in Sec. IVI| . Further- 
more, theoretical and experimental NMR studies on 
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fullerenea^ have been carried out as well as ab initio 
calculations on hyperfine interaction in small graphene 
flakes^ An analytical investigation of nuclear-spin in- 
teractions of electrons confined to graphene and CNTs, 
however, has so far been missing. 

In our work, we analytically calculate the interaction of 
a single electron confined to a graphene ribbon or a CNT 
(see Fig. [1} with the spins of the surrounding 13 C nu- 
clei. For CNTs we find an interesting interplay between 
isotropic and anisotropic hyperfine interactions which de- 
pends on the CNT geometry and which leads to a highly 
anisotropic Knight shift and an unusual alignment of the 
nuclear spins around a CNT circumference (in the ground 
state). Furthermore, we calculate the decoherence dy- 
namics of the electron spin and find that, even without 
manipulating the nuclear spins, the hyperfine-associated 
decoherence times can be on the order of tens of microsec- 
onds or longer, depending on the relative abundance of 
13 C nuclei. These timescales are much longer than the 
ones typically found for III-V semiconductor quantum 
dots, making carbon-based quantum dots promising spin- 
qubit candidates. 

Our paper is organized as follows: In Sec. |TT] we de- 
termine the wavefunction of a conduction-band electron 
in a CNT. In Sec. IIIII we derive an effective Hamilto- 
nian for the nuclear-spin interactions and find a highly 
anisotropic coupling of the electron to the 13 C nuclear 
spins, which results in an anisotropic Knight shift and 
an unusual alignment of the nuclear spins around the 
nanotube, as shown in Sec. IIVI In Sec. [V] we deter- 
mine the electron-spin dynamics and the associated spin 
decoherence times. We give a detailed comparison of our 
results with previous studies on comparable systems in 
Sec. |VH and conclude in Sec. IVIII 

II. BONDS AND BANDS 

The carbon atom features six electrons: two core elec- 
trons and four valence electrons. In a solid, the four 
valence electrons (which are in 2s and 2p states) form 
the bonds with the nearest-neighbor atoms. In two- 
dimensional graphite, or graphene, each atom has three 
nearest neighbors, and three of the four valence electrons 
form sp 2 -hybridized bonds (called tr-bonds) with those 
neighbors, while the fourth electron is in a so-called w- 
state perpendicular to the cr-bonds^ The 7r-electrons in 
graphene determine the band structure near the Fermi 
energy, while the cr-electrons form more remote bands. 
A conduction-band electron is therefore in a 7r-state. 

A CNT can be regarded as a graphene sheet that had 
been rolled up along some direction c, defining a symme- 
try axis t of the CNT. The circumferential and transla- 
tional vectors are defined in terms of the basis vectors ai 
and a 2 as (see Fig. [2]) 

c = n&i + ma 2 , t = t\B.i + £2 a 2- (1) 

Here, n,m G No are the chiral indices, and t\ = (2m + 




FIG. 2: Definitions used throughout this work: ai and &2 
denote the basis vectors of the honeycomb lattice, Rj are the 
relative nearest-neighbor positions, c and t are the circum- 
ferential and transverse vectors, respectively, and 6 = 9 nm is 
the chiral angle between c and ai. 

n)/du, t% — —(2n + m)/d,R with cIr being the greatest 
common divisor of 2n + m and 2m + n. We denote the 
chiral angle between c and ai by 9 = 9 nm . 

The curvature of the CNT causes a geometrical tilt- 
ing of the cr-bonds which results in an sp-hybridization 
of the conduction-band states. In Ref. [39|, this effect has 
been studied perturbatively in lowest order in the small 
parameter 2tt/L, where L = \/n 2 + nm + m 2 is the cir- 
cumference of the CNT in units of the lattice constant 
a -2.5 AM. 

The wavefunction of a conduction-band electron in a 
periodic crystal is given by Bloch's theorem: = 
^=e lk r Ukcr(r), where Na is the number of atomic sites 
in the crystal and the Bloch amplitude itko-(r) has the 
periodicity of the lattice. Following Ref. |4fJ, we ap- 
proximate the Bloch amplitude u^ a (v) at the K and K' 
points by a linear combination of hydrogenic orbitals, 
Ua{r) = J2r 7r ( r — R)j where 3 ^ 

7r(r) =N nm {ip 2 pA r ) + — ^-^ 2s ( r ) 

+ frfn(3CJ V>2 Pt (r) + cos(3CJV>2p c (r)) } , (2) 

and the sum runs over all lattice sites in the CNT. In 
the above, -02s represents a hydrogenic 2s orbital, ip2p t , 
ip2p a and ?/>2pj_ are, respectively, the 2p orbitals along 
the transverse, circumferential, and radial direction of 
the CNT, and N nm normalizes the Bloch amplitude to 
two atoms per unit cell. 6' nm is the angle between c 
and Ri, and we can write sin(3#^ m ) = (n — m)(2n 2 + 
5nm + 2m 2 ) /2L 3 and cos(36>; m ) = 3V3nm(n + m)/2L 3 
in terms of the chiral indices.— The hydrogenic orbitals 
consist of a radial and an angular part, e.g., ip2s( r ) — 
i?2o( r )^o°(^; V 5 )- We will choose a local coordinate sys- 
tem at each lattice site, such that i/)2p t , ^2p a , and ?p2p ± 
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correspond to hydrogenic 2p orbitals in z-, y-, and x- 
direction, respectively (see below). The radial compo- 
nents of the hydrogenic orbitals depend on an effective 
screened nuclear charge Z B g 'seen' by the electron^ 

Our choice of the Bloch amplitude implicitly assumes 
that the electron is tightly bound to the nuclei, i.e., that 
the radial component of 7r(r) drops off fast on the scale of 
the nearest-neighbor distance. We have estimated that 
\ir(r + R n n )| 2 /|7r(r)| 2 ~ 10 -3 for any nearest-neighbor 
lattice vector R n .n., justifying our assumption. 

In a quantum dot, the electron is delocalized over many 
lattice sites and its Bloch amplitude is modulated by an 
envelope function <& CT defined by the confinement poten- 
tial (see Ref. Hi for the envelope function of a quan- 
tum dot defined by a rectangular confinement poten- 
tial in a semiconducting CNT). Including the spin states 
\ a ) = I T> 4)j we write the electron states as 

|f CT ) = |$ CT ;n CT )|a), (3) 

where, in the envelope-function approximation, 
(r\$ a ;ucr) = <f> cr (r)w cr (r). Since we did not include 
spin-orbit interactions in our model, the orbital wave- 
function is independent of the spin state, and we may 
omit the subscript a: $ CT (r) = $(r) and u a (r) = it(r). 



III. NUCLEAR-SPIN INTERACTIONS 
A. Carbon nanotubes 

There are three terms that couple the spin of the con- 
fined electron to the nuclear spins in the CNT: the Fermi 
contact interaction, the anisotropic hyperfine interaction 
and the coupling of electron orbital angular momentum 
to the nuclear spins. These interactions are represented 
by the Hamiltonians^ 2 - 
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respectively, where 7s = 2fi B , Jj k = 9 ]k HN, Ps is the 
Bohr magneton, gj k is the nuclear g-factor of isotopic 
species j k , Pjv is the nuclear magneton, p is the vacuum 
permeability, = r — is the electron-spin position 
operator relative to the nucleus, d ~ Z x 1.5 x I0~ 15 m 
is a length of nuclear dimensions, Z is the charge of the 
nucleus, and n k — v k jr k . S and Lfe = r k x p denote 
the spin and orbital angular-momentum operators (with 
respect to the fc th nucleus) of the electron, respectively. 
The cutoff 1 + d/ru comes from the Dirac equation (see, 
e.g., Ref. |42| ) and avoids unphysical divergences from 
expectation values of the Hamiltonians h\ and /ig. In 
the problem considered here, this cutoff may be omitted 



for the following reasons: (i) The expectation values of h\ 
and with respect to an s-state vanish identically due 
to the spherical symmetry of the wavefunction and the 
vanishing orbital angular momentum, respectively. The 
expectation values with respect to a p-state are non-zero, 
but the p-wavefunction goes to zero sufficiently fast at 
the position of each nucleus, thus avoiding a divergence, 
(ii) As mentioned above, the electron wavefunction does 
not extend significantly to the nearest-neighbor lattice 
sites. Hence, within our tight-binding approximation, the 
orbital 7r(r) centered around some nucleus cannot cause 
a divergence at the position of a nearest neighbor. 

We note that the orbital angular momentum which 
appears in Eq. ^ is associated with the electronic mo- 
tion around the k th nucleus and is described by the Bloch 
part of the electron wavefunction. On the other hand, it 
has been shown in Ref. |43| . that interatomic currents can 
occur along the nanotube circumference to which another 
orbital angular momentum C may be associated which is 
described by the envelope part of the electron wavefunc- 
tion. However, the hyperfine coupling strength is defined 
via the Bloch part of the electron wavefunction and there- 
fore a consideration of the coupling between the valley 
degrees of freedom via envelope-function-associated an- 
gular momenta is beyond the scope of the present work. 

We will consider CNTs and graphene with different 
abundances of the nuclear isotopes 12 C and 13 C. While 
12 C does not carry a nuclear spin, the nuclear gyro- 
magnetic ratio of 13 C is non- vanishing and given by 
7 i3 C = 7.1 x KT 27 J/T. 

The Fermi contact interaction (0| yields a finite con- 
tribution for s-states, but vanishes for p-states. The 
anisotropic hyperfine interaction ([5]) and the coupling 
of orbital angular momentum ^ vanish for s-states be- 
cause of their spherical symmetry and zero orbital an- 
gular momentum, but yield a finite contribution for p- 
states. Therefore, when considering CNTs, all three in- 
teractions (JU) - © have to be taken into account be- 
cause of the sp- hybridized electron states ©, while for 
graphene, only the interactions ([5]) and (J6j) are relevant, 
due to the purely p-type wavefunction (corresponding to 
the limit n,m — > 00 in Eq. ([2])). 

We first calculate matrix elements of the interactions 
(HJ) - © with respect to the electron wavefunction ©, 
which will lead to effective spin Hamiltonians and to the 
associated coupling strengths in the CNT case. Through- 
out this section, we will consider a CNT that consists 
only of spin-carrying 13 C isotopes. The possibility of dif- 
ferent nuclear isotope abundances will then be taken into 
account in Sec. [V] From the CNT results, it will be pos- 
sible to perform the 'graphene limit', which we postpone 
to Sec. IHTBI 

We start with the Fermi contact interaction and cal- 
culate 
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assuming that the electron-spin density does not depend 
on the lattice site, which is justified if, e. g., the envelope 
function $(r) describes the ground state of the quantum 
dot. The effects of a site dependence of the electron-spin 
density have been recently considered in Ref. |44j. Eval- 
uating the spin matrix elements leads to the following 
effective spin Hamiltonian: 



s • I, 



(8) 





(10, 0) 


(20,0) 


(5, 10) 


(9, 15) 


(5,5) 


(10,10) 


(20, 20) 


A 1 


0.19 


0.05 


0.12 


0.05 


0.26 


0.07 


0.02 


A% 


0.59 


0.60 


0.59 


0.60 


0.58 


0.60 


0.61 


Al 


-0.29 


-0.30 


-0.30 


-0.30 


-0.29 


-0.30 


-0.30 


A% 


-0.30 


-0.30 


-0.30 


-0.30 


-0.30 


-0.30 


-0.30 



TABLE I: Estimated hyperfine coupling strengths (in /xeV) 
for nanotubes of different chiralities (n, m). 



with coupling constants A*jp = AiVo\$(rk)\ 2 (where vq 
is the volume of a primitive unit cell) and the associated 
coupling strength 



Al = 3^1 Ar "" 



A. 



2 

nm ' 



(9) 



where we have introduced (3 nm = w/2^/SL, and do is 
the Bohr radius. The normalization factor N nrn can be 
determined by normalizing Eq. ([2]) to two atoms per unit 
cell: 



N — 2 

1 v n:m, — ^ 



L' 2 



4L 2 



(10) 



We have evaluated Eq. §§§ for CNTs of different chirali- 
ties in Table [U For typical CNTs, the coupling strength 
A\ is about three orders of magnitude smaller than that 
for an electron in a GaAs quantum dot (A^ aAs ~ 90/ieV, 
see Ref. l45h . for two reasons: (i) The hybridization pref- 
actor N nm f3 nm is on the order of 0.05 and enters quadrat- 
ically into A\. (ii) The effective nuclear charge (which 
enters in third power into A\) is Z^ s ~ 3.2 for carbon, 
but Zfi ~ 7.1 for gallium and ~ 8.9 for arsenical 

The sign of the isotropic interaction is positive for 
all nanotube diameters, in contrast to results reported 
previously^ where spin relaxation of conduction elec- 
trons has been calculated numerically and a sign change 
of the hyperfine coupling constant A\ has been found for 
small nanotube diameters. This is due to the fact that in 
Ref. [46| the coupling of the Is core electrons to the nu- 
clear spins is included, which, however, is irrelevant for 
the hyperfine interaction of a conduction electron and, 
hence, for the electron spin dephasing considered in this 
article. 

Now we look at the anisotropic hyperfine interaction. 
Due to symmetry, the s-part of the hybridized wave- 
function does not contribute. Taking matrix elements 
(J5>„ 1 h 2 | ^! a > ) j us t like above, we arrive at an effective 
Hamiltonian 



k 

,(2j) 



(11) 



with coupling constants A},' 3 ' = A 2 vo\Q(rk)\ 2 and the 
coupling strengths 



2 " 120vrag nm j ' 



(12) 



where 



A* = 1 



2 ' nm i 



(13) 



\ = -\- \$L sin 2 (3C„) + 0tm cos 2 (3e' nm ), (14) 

A 2 = -\ + Plm sin 2 (3C„) - \ film cos 2 (3^ m ). (15) 

We see that the coupling induced by the anisotropic hy- 
perfine interaction is different in all three spatial direc- 
tions. Recall that we have labeled our axes such that z, 
y, and x refer to the translational, circumferential, and 
radial directions, respectively. We show typical values for 
A 2 in Table |U The direction of strongest hyperfine inter- 
action is radial to the nanotube^ It is interesting to note 
the competing signs: the Fermi contact interaction (ex- 
pressed via A\) enhances the anisotropic hyperfine inter- 
action along the radial direction, but reduces it along the 
circumferential and transverse directions. Furthermore, 
the signs of the hyperfine coupling along different direc- 
tions indicate ferro- and antifcrromagnetic alignment of 
the nuclear spins with respect to the electron spin (in the 
ground state). We will come back to this in Sec. IIVI 

Finally, we address the coupling of electron orbital an- 
gular momentum. Calculating matrix elements of it is 
straightforward to see that this interaction vanishes iden- 
tically by applying the operators Lj to the hydrogenic or- 
bitals appearing in Eq. ((2|), then evaluating a \h^\^ a >) 
and summing up all contributions. 



B. Graphene 

We consider the 'graphene limit' corresponding to let 
n, m — > oo in all expressions in Sec. IIII Al Then 
Pnm — ► and N nm — > 1 and, denoting quantities related 
to graphene with a tilde, 



Ai = 0, Al 



A', 



240vra^ 



(16) 



Inserting numbers, this yields A 2 = A 2 ~ —0.3 fieV and 
A 2 ~ 0.6 fieV. Note that in our notation the ^-direction 
is perpendicular to the graphene plane. 

Surprisingly, these numbers are not much different 
from those estimated for A 2 in Sec. IIII Al Naively, one 
might have expected a weaker hyperfine interaction in 
graphene (as compared to CNTs) due to its flatness and 
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FIG. 3: Alignment of the nuclear spins (in the ground state) 
due to the anisotropic Knight shift. The electron is assumed 
to be prepared in the eigenstate of 5™ with eigenvalue +1/2, 
see Eq. (|19[1 . and points along x at each nuclear-spin site. 



are small with respect to the energies associated with the 
Knight shift: k B T < (A 1 + A J 2 )/N ~ lpeV...lneV, 
where ks is the Boltzmann constant. We write the 
Hamiltonian describing the hypcrfine interaction of one 
nucleus with the delocalized electron as 

H' M = S T AI = S T Al, (17) 

where S T = (S x , 5", S z ), A = di&g(A x ,A y ,A z ), I = 
(I x ,iy,I z ) T , S T = S T R\ A = RAR\ and I = RI, 
with the rotation R given by 

/ cos C — sin C 0\ 
R = sin( cos( . (18) 



The operators carrying a tilde hence describe the inter- 
action in the global coordinate system (x,y, z), and the 
coupling tensor is given by 



the vanishing contact interaction ([JJ. As it turns out, 
however, even for small CNTs the contact interaction is 
only a small correction to the anisotropic hyperfine in- 
teraction ([5]), so that the latter is the main hypcrfine 
contribution for both CNTs and flat graphene. 



IV. HYPERFINE-INDUCED ANISOTROPIC 
KNIGHT SHIFT 

In Sees. |TT] and IIIH we have formulated the hypcrfine 
problem in terms of a local coordinate system at each nu- 
cleus. In this section, we want to look at the Knight shift 
of the 13 C nuclear spins due to hyperfine interaction with 
the conduction electron. The isotropic Knight shift due 
to interaction with both the sp-hybridized conduction- 
band electron and the Is core electron has been studied 
in Ref. H?]. Electron-spin resonance spectra of 13 C in 
7r-electron radicals have been analyzed in Ref. |48| . 

From our considerations in Sec. IIII Al it is clear that 
the Knight shift induced by the anisotropic hyperfine in- 
teraction (~ A J 2 /N) exceeds the isotropic Knight shift 
(~ Ai/N) by roughly one order of magnitude (N is the 
number of nuclei in the dot). The Knight shift in CNTs 
(and graphene) is therefore highly anisotropic. 

We introduce the following global coordinate system: 
x — xcosC — ysinC, y = xsinC + ycos£, z = z, such 
that the iy-plane cuts out a cross section of the CNT 
and C is the coordinate describing the position on this 
cross section (see Fig. [3]). 

The electron is assumed to be prepared in a fixed state 
and to be delocalized over the CNT cross section. At each 
lattice point occupied by a 13 C nucleus, the nuclear spin 
will align itself in such a way that the hyperfine energy is 
minimized in the ground state of the nuclear spins. This 
ground state can only be achieved for temperatures that 



(Ax cos 2 C + A y sin 2 C (Ax - A y ) sin C cos C \ 
(A x - Ay) sin C cos ( A x cos 2 ( + A y sin 2 CO. 
A Z J 

Here, Aj = Ai+A J 2 is the sum of isotropic and anisotropic 
hyperfine couplings along the (local) j-direction (j — 
x,y,z), see Eqs. (|9|) and (fT2)) . Recall that A x > and 
A y ,A z < for CNTs (see Sec. IIII Aft . For an electron 
spin in an eigenstate of S x , the interaction reads 

//;;, ■ A. r: ,r<iS> (19) 

with Axx = A x cos 2 C + A y sin 2 C and A X y — (A x — 
A y ) sinC cosC- 

In Fig. [3] we show the alignment of the nuclear spins 
(assuming the nuclear spins to be in their ground state) 
due to the anisotropic Knight shift induced by the hyper- 
fine interaction with a conduction electron whose spin is 
prepared in the eigenstate of S x with eigenvalue +1/2, 
i.e., pointing along the x-direction at each nuclear site. 
We have assumed that the electron is evenly distributed 
around the CNT cross section, which can be seen to be 
justified from the envelope functions calculated in Ref. [ID 
for semiconducting CNTs subject to a rectangular con- 
finement potential. We observe an interesting interplay 
between ferro- and antiferromagnetic coupling along the 
two spatial directions, which is a direct consequence of 
the CNT geometry and the strong anisotropy of the hy- 
perfine interaction. In particular, the hyperfine interac- 
tion does not vanish when we average over the CNT cir- 
cumference. This could lead to other interesting effects. 
For instance, it has been shown in Ref. that in a Lut- 
tinger liquid, a non-vanishing average hyperfine field can 
lead to a transition into a helically ordered phase (along 
the tube axis) of the nuclear spins below some critical 
temperature. 
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V. ELECTRON-SPIN DECOHERENCE 
A. Carbon nanotubes 

Based on our analysis in Sec. IIIII the electron-spin dy- 
namics in a CNT arc described by the following hyperfine 
Hamiltonian: 

ff hf = h-s, h, Y^Ali:. (20) 

k 

where A\ = A,-i; |<I>(r fe )| 2 and A } — A 1 + A\, see Eqs. 
© and fl2]). 

We assume that an external magnetic field B z has 
been applied along the symmetry axis of the CNT (see 
Fig. Q] (b)), and that the induced Zeeman splitting be- 
tween the spin states is larger than the hyperfine coupling 
strength: b = gfisBz > Aj. Assuming g ~ 2, this cor- 
responds to very moderate fields of B z > 5 mT. Within 
this limit, relaxation-induced decoherence is suppressed 
by the small parameter Aj/b, and the main source of 
decoherence is pure dephasing due to nuclear-field fluc- 
tuations along the CNT symmetry axis. The relevant 
Hamiltonian is given by the z-part of Eq. ([20| and the 
Zeeman term: 

H = {b + h z )S z . (21) 

Assuming that all nuclei carry non-zero spin, the dynam- 
ics of the transverse spin is then given by a Gaussian de- 
cay with timescale given by y/W /A z , where N is the total 
number of nuclear spins in the quantum dotJ^ 

One advantage of carbon-based nanostructures is the 
low natural abundance of the spin-carrying isotope 13 C. 
In order to investigate this advantage, we allow for ar- 
bitrary 13 C abundances: we denote the total number 
of nuclei in the quantum dot (defined via the enve- 
lope function, see below) by N, and by AT 12 and N i3 
the number of 12 C and 13 C nuclei, respectively, such 
that N — N\2 + Nis. This generalization has two rel- 
evant effects: (i) Summing over all N nuclei, we get 
Sfc A? ~ V As j i- e -j the total hyperfine coupling is weak- 
ened by the ratio r] = N13/N. (ii) The polarization p of 
the nuclear-spin system (0 < p < 1) is determined only 
by the distribution of spin- up and spin-down 13 C nuclei, 
while being unaffected by the 12 C nuclei. 

For ./V13 > 1 we can use the central limit theorem 
(compare with Ref. [Toh which yields the following Gaus- 
sian dynamics for the transverse electron spin (written 
in the frame rotating with frequency (ui + pr]A z /2)/h, 
where uj = b — b^ with the nuclear Zeeman energy 
bN = 9nHnB z ): 

(S+) t = (S + )oe- t2 /< (22) 

Here, S + = S x + iS v , with the electron-spin operators 
5 3 = (Tj/2 and the Pauli operators aj. The characteristic 
timescale for the decay is given by^ 

' y/T^p 1 Vn^a z 




I 1 , , , 1 , , , 1_ 

0.1 0.5 0.9 

relative 13 C abundance 



FIG. 4: Electron-spin decoherence time r c in a (20,0) CNT as 
a function of the relative 13 C abundance N13/N, under the 
condition of a magnetic field B z > 5 mT along the symmetry 
axis of the CNT. We have chosen a completely unpolarized 
nuclear system (p — 0) and iV = 6 x 10 5 , see text below 
Eq. J23}. Inset: Electron-spin dynamics in CNTs with 13 C 
abundances of 1% (solid blue curve) and 99% (dashed red 
curve), in the frame rotating with u + pNis, see Eq. (1221) . 



We can see from Eq. (|23|) that the decoherence time r c 
has an interesting non-linear dependence on the abun- 
dance of 13 C nuclei. 

The total number N of nuclei in the quantum dot is 
determined by the quantum-dot confinement. Assuming 
a rectangular confinement along the symmetry axisi^ and 
a (20,0) CNT of length 300 nm, we estimate N ~ 6 x 
10 5 . We show the electron-spin decoherence time t c as 
a function of the relative 13 C abundance in Fig. 2] for a 
completely unpolarized nuclear bath (p = 0). We see that 
decoherence times of several tens of microseconds can be 
expected for higher 13 C abundances. For the natural 13 C 
abundance of about 1%, we estimate r c > 200 /is from 
Eq. ([23]) . The inset of Fig. [4] shows the Gaussian decay 
of the spin coherence (Eq. P?2"]) for CNTs containing 1% 
(solid line) and 99% (dashed line) 13 C nuclei. We note 
that further increase of the decoherence time r c , say by 
a factor of x, would require isotopic purification and the 
reduction of the natural 13 C abundance by about a factor 
of x 2 . The decoherence law in Eqs. and P3| of 

course breaks down at the point when only a few nuclear 
spins are present in the dot. 



B. Graphene 

For an electron confined to a graphene quantum dot, 
the contact interaction vanishes identically {A\ = 0, see 
Eq. (|16ph and only the anisotropic hyperfine interac- 
tion contributes to spin decoherence. We assume an ex- 
ternally applied magnetic field B x perpendicular to the 
graphene plane (see Fig. H (a)), giving rise to a Zeeman 
splitting b — g±fisB x , where g± ~ 2 is the electron g- 
factor in the out-of-plane direction. If the Zeeman split- 



7 



'Si 



40 



20 




0.1 0.5 0.9 

relative 13 C abundance 

FIG. 5: Electron-spin decoherence time t c in graphene as 
a function of the relative 13 C abundance N13/N, under the 
condition of a magnetic field B x > 5 mT perpendicular to 
the graphene plane. We have chosen p = and N = 4 x 
10 5 (see text in Sec. IVB|I . Inset: Electron-spin dynamics in 
graphene with 13 C abundances of 1% (solid blue curve) and 
99% (dashed red curve), in the rotating frame. 



ting is much larger than the energy associated with the 
transverse hyperfine terms, b 3> A%, A\, the electron-spin 
dynamics are governed by a Gaussian decay similar to the 
CNT case (|22|) . but with a characteristic timescale given 
by 



2/t 



N 



V 7 ! 



(24) 



For a quantum dot with width FK = 30 nm and length 
L = 30 nm, we estimate the total number of nuclei to be 
N ~ 4 x 10 5 . We show f c as a function of the relative 
13 C abundance in Fig. [5j for higher 13 C abundances, 
decoherence times f c > 1/xs can be expected, while for 
the natural 13 C abundance of 1%, we estimate f c > 80/xs. 



C. Comparison and discussion 

It is interesting to note that the electron-spin deco- 
herence times in graphene are shorter than in CNTs. 
Naively, one would assume that the CNT curvature and 
the associated hybridization would lead to an enhance- 
ment of the nuclear-spin interactions due to the contact 
interaction As it turns out, however, the contact in- 
teraction in a CNT has a competing sign as compared 
to the anisotropic hyperfine interaction along the CNT 
symmetry axis and, hence, effectively reduces the total 
hyperfine coupling strength along the direction of the ex- 
ternal magnetic field. For ultra-small CNTs it might hap- 
pen that the total hyperfine coupling along the symmetry 
axis approaches zero^ This, however, is beyond the va- 
lidity of our theory. 

For typical CNTs, the dominant contribution to 
the nuclear-spin interactions in CNTs comes from the 



anisotropic hyperfine interaction. This is because the 
amount of s-orbital admixture in the hybridized wave- 
function of a conduction-band electron is rather small: 
on the order of a few percent. 

In our considerations throughout this section, we have 
neglected the hyperfine terms which are transverse to the 
externally applied magnetic field, i.e., radial and circum- 
ferential in the case of a CNT, and in-plane in the case 
of graphene. This gives a good first approximation of the 
decoherence time, as long as the external magnetic field 
is large enough to suppress spin flips induced by the hy- 
perfine terms transverse to the external field. For both 
CNTs and graphene, magnetic fields of order B > 5 mT 
are sufficiently strong to achieve this. 



VI. COMPARISON WITH PREVIOUS WORK 

In Table [TT] we compare our results for the hyperfine 
interaction in graphene and CNTs with those given in 
earlier publications for comparable systems. Yazyev's 
values22 for the hyperfine interactions in small graphene 
flakes were derived from DFT calculations including the 
(T-bands. The values for the anisotropic hyperfine inter- 
action are similar to our results but a non- vanishing (and 
even negative) value for the isotropic hyperfine interac- 
tion is reported, which is, however, associated with the 
coupling between the nuclear spins and the Is core elec- 
trons, which is irrelevant for the dephasing of the conduc- 
tion electron. Pennington and Stenger— give estimates 
for the hyperfine interaction in fullerenes and report cou- 
pling constants which are slightly larger than our values. 
The relatively large value for the isotropic hyperfine in- 
teraction can be explained by the stronger curvature in 
Ceo molecules, leading to stronger sp-hybridization of the 
electron states, as compared to CNTs. Goze-Bac et al^ 
have estimated the hyperfine interaction in CNTs based 
on measurements of the chemical shift of 13 C. 

Churchill et aL— have estimated the hyperfine in- 
teraction in carbon nanotubes from transport measure- 
ments in double quantum dots, but do not comment on 
the anisotropy of the interaction. A hyperfine coupling 
strength of ~ 100 peV for pure 13 C CNTs is reported 
in their work, in clear contrast to our results. Cur- 
rently, this discrepancy is not understood. One possi- 
ble explanation^ 1 - for this might be that the theory de- 
veloped in Ref. Km for standard GaAs quantum dots 
and used in Ref. [34j to deduce the hyperfine coupling 
strength gets modified by the valley degeneracy occurring 
in CNT and thus this could lead to different conclusions. 
Another possible explanation 34 could be that the effec- 
tive electron-nuclear interaction gets greatly enhanced 
by electron-electron interactions and the one-dimensional 
character of the system. A similar renormalization was 
recently noticed in the context of nuclear magnetic order- 
ing in CNT42 Clearly, this is an interesting open prob- 
lem, which, however, requires separate study. 
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Ai [fJ,eV] 


A% [/LteV] 


* 11 Z r t rl 

AY \peV] 


our values (CNT) 


0.05 


0.6 


-0.3 


our values (Graphene) 





0.6 


—0.3 


Yazyev 3 ' (Graphene flakes) 


-0.2 


0.6 


-0.3 


Pennington, Stenger~ (Ceo) 


0.1 


0.9 


-0.5 


Goze-Bac est al^ (CNT) 


0.04 


0.9 


-0.5 



TABLE II: Comparison of our hyperfine coupling strengths 
with the values from previous publications for comparable 
systems. We give our values for a (20,0) zigzag nanotube and 
for graphene. See Sec. I VII for a detailed discussion. 



VII. CONCLUSIONS 

We have calculated the nuclear-spin interactions and 
the resulting spin dynamics of an electron confined to 
a CNT or graphene quantum dot. In graphene, only 
the anisotropic hyperfine interaction couples the elec- 
tron to the nuclear spins, due to the purely p-type elec- 
tron wavefunction. In a CNT, curvature induces an sp- 
hybridization of the electron orbital, opening a new chan- 
nel of spin decoherence via the Fermi contact interaction. 
However, for typical CNTs, the Fermi contact interaction 
is only a small correction to the anisotropic hyperfine in- 
teraction, the latter being the main source of nuclear- 
spin-induced decoherence of the electron spin. We found 
the total hyperfine coupling strength of an electron with 
the 13 C nuclei to be less than 1 /LteV for both graphene 
and CNTs quantum dots - about two orders of magni- 
tude smaller than the hyperfine interaction of an electron 
in a GaAs or In As quantum dot. 

We have used a simple model for the sp-hybridization 
in CNTs, from which we have derived the hyperfine in- 
teraction. We have checked, however, that a numerical 
tight-binding bandstructure calculation yields hybridiza- 
tion on the same order of magnitude as the geometrical 
approach used in this worki^ 3 - Nevertheless, a more rigor- 
ous bandstructure calculation including also the influence 
of the CNT CT-bands would be desirable in the future. 

For CNTs, we have found an interesting interplay 
of hyperfine couplings along different spatial directions, 
leading to a highly anisotropic Knight shift and an align- 
ment of the 13 C nuclear spins around the CNT in the 
ground state of the nuclear spins. This result is particu- 
larly interesting when viewed in context of the hyperfine- 
induced nuclear phase transition predicted in Ref. for 
a Luttinger liquid, which requires a non- vanishing mean 
value of the hyperfine field. The strong anisotropy of the 



hyperfine field that we have found in our work thus gives 
further evidence that such phase transitions could occur 
in CNTs. We emphasize that this anisotropy is present 
in CNTs of any chirality, in particular also in metallic 
CNTs. The results presented here have been used very 
recently in the context of transport though CNT quan- 
tum dots4^ 

Furthermore, we have estimated typical electron-spin 
decoherence times in CNT and graphene quantum dots. 
We have shown that relaxation-induced decoherence due 
to nuclear spins can be suppressed by applying a mag- 
netic field of order 5 mT to the system, leaving only pure 
dephasing due to fluctuations of the nuclear magnetic 
field. We have estimated that for a (20,0) zigzag CNT 
quantum dot containing N ~ 6 x 10 5 nuclei, and for a 
magnetic field applied along the CNT symmetry axis, the 
associated decoherence time is of order r c > 1 fis, depend- 
ing on the relative 13 C abundance. We emphasize that 
our analytical treatment of the hyperfine problem applies 
to CNTs of any chirality. For a graphene quantum dot 
containing JV = 4 x 10 5 nuclei, and for a magnetic field 
applied perpendicular to the graphene plane, the deco- 
herence time is of order f c > 0.5 /us, again depending on 
the relative 13 C abundance. 

The hyperfine interaction in the systems we have con- 
sidered here is rather weak. Therefore, it could, in prin- 
ciple, be that other mechanisms, such as spin-orbit inter- 
actions, limit the lifetime of spin-state superpositions on 
timescales comparable to those we have estimated here 
(see Ref. [H for details) . 

The decoherence times we have estimated throughout 
this work are among the longest reported so far. This 
makes quantum dots based on carbon materials attrac- 
tive spin-qubit candidates. In particular, the tunability 
of the average hyperfine coupling strength via the abun- 
dance of spin-carrying 13 C nuclei could be used to achieve 
an optimal balance between a long electron-spin decoher- 
ence time and a sufficiently strong coupling to control the 
electron-spin state by manipulating the nuclear-spin sys- 
tem. 
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